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Abstract 


I'his report presents the prediction of zenith runge refraction from surf an 
measurement* of meteorological parameters. R»*fr activity is separated into wet 
(water vapor pressure) ami dry (atmospheric pressure) components The Integra 
tion of dry refractivity is sliown to be exact. Attempts to integrate wet refractivit> 
directly prove ineffective; however, several empirical models developed l>y the 
author and other researchers at Jl’l. are discussed The best current wet refraction 
model is here considered to be a separate day/night model ("Berman (l)/N 
which is proportional to surface water vapor pressure and Inversely proportional 
to surface temperature. The standard deviation of this model is considered to he: 

a ~ 1.5 — 2.0 cm 

Methods are suggested that might improve the accuracy of the wet range retrac- 
tion model, however, the information content in surface parameters is considered 
insufficient to allow a surface measurements model for wet range refraction to 
result in a standard deviation lower than 

1.0 1.5 cm 
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The Prediction of Zenith Range Refraction From Surface 
Measurements of Meteorological Parameters 


I. Introduction 

In the last two decades, increasingly sophisticated deep 
space missions haw placed correspondingly stringent re- 
quirements on navigational accuracy. As part of the effort 
to increase navigational accuracy, and lienee the quality 
of radiometric data, much effort has been extended in an 
attempt to understand and compute the tropospheric 
effect oi. range (and hence range rate) data The general 
approach adopted has Ih-cii that of computing a zenith 
range refraction, and then mapping this refraction to any 
arbitrary elevation angle via an empirically derived func- 
tion of elcsutioii. Thus if 


The relevant parameters necexsar) to the determination 
of tro|K»pheric range refraction are as follows: 

n index of refraction 

V rcfractivity 

z height al>o\’c station, km 

r range, km 

f refracted range, km 

r speed of light, km/s 

I time, s 


sR 

zenith range refraction, cm 

where 

n 

elevation angle, deg 

n = 1 + 10 " X 

f(») 

elevation “mapping” function 

One begins b> considering the physical signal path range 

then 


as the troposphere is traversed: 

ilz 

SR. 

sR {/(#)) 

(It — 

c 

where 


and 

/(90) = 1 


r — fedt fdz 
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X 


t 

I 


I 

[' 


One then considers (in apparent refracted range f ("sig- 
nal retardation") with signal velocity i> (c > o): 


where 

ND = "dry” rcfractivity 


o 




t It' 



and 


N\V = "wet” rcfractivity 


C,C, 



r'^J cdt* c ~ Jn dz 
so that the corresponding range refraction would be 


and 


C, -= 0.770 
C, =-1810.0 


ARra 10* [r'- r] 

•- 10 1 [/n dz — fdz] 

= 10’/(» - 1} dz 
= 10“/10-Wr/s 
= lO-'/Wrfs 

Prior to 1070, attempts to calculate this quantity generally 
assumed that 

N(z) = N, exp(-B,s) 

where 

= surface rcfractivity 
B t = “inverse scale height,” km (si 0,1) 


P = (totnl) atmospheric pressure, N/m s {1 mbar 

= 10* N/m 3 } 

T = temperature, K 

PW — water vapor pressure, N/m 3 

/ AT — 11 \ 

= 010 {n/7} exp y- f] , _ c j 

IUI — relative humidity (1.0 = 100%) 

A = 7,4475 In (10) = 17.1-185 
B = 2034,28 In (10) = -1684.1 
C = 38.45 

At that time it was shown that the quantity; 


with Z3 t an empirically determined constant, The total 
range refraction (considering that tropospheric type range 
refraction is nil by about 75 km) was then 



could be integrated exactly, while the quantity 


AR = 10- 1 / Ng exp (—Biz) dz 


could be integrated approximately, which represented a 
considerable breakthrough since the dry component of 
rcfractivity contributes the vast bulk of the refractive 
effect, viz, 


^ ' jo ( " B > 3 ) d z 


10-> 


N. 

lit 


/: 


NW(z)dz 


In 1070, this author broke new ground (Ref. 1) by con- 
sidering the wet and dry components of rcfractivity sepa- 
rately, viz, 


N-ND + NW 


0.00 3 


f* ND(z) dz 

J” N(z) dz 


30.99 
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I 




I 


and 


where 


I NW(x) dz 

0.01 Z-ljz — £ 0.10 

/ N(z),U 

Thi* rr|Mirt rcviewi tin work leading to the exact inte- 
gration of the dry component ol refractivity (Section III, 
and then present* a theoretical integration ol the wvt 
eoni|xinent of refractivity. a* well a* more recent work hy 
the author and other researchers in tin* field towards the 
refinement of an expression for wet refraction (Section 
III). Finally. Section IV presents a summary, and com- 
ment* on future prospects for additional gain* in accuracy 
in tin* determination of range refraction from the mea- 
surement of surface meteorological parameters. 


II. Determination of Dry Zenith 
Range Refraction 

Dry range refraction is given hy the following 
expression: 



SD(z)dz 


where 


Afti = “dry" range refraction, cm 

In tin* following Subsection*, it will In* shown that 
different atmospheric assumptions all h ad to tin* same 
results for the determination of dry range refraction. 

A. Isothermal Atmosphere 

As was mentioned in Section I. early attempts to inte- 
grate (total) refractivity assumed a refractivity profile of 

.V(z) = N.exp(-fl,z) 

If one limits the discussion to dry refractivity only, one 
finds that this assumption is identical to the assumption of 
an isothermal, or constant temperature, atmosphere. The 
definition of an isothermal atmosphere (sec Ref. 2, p. 82) 
is as follows: 


T. surface temperature. K 

P. surface pressure, N'/m* 

U gravitational acceleration. USOfl < m/» 
H IYru*ct C .a* Constant <>2S7 J/gk 

~ 94.1 K/km 


Calculating dry rang** refraction one* would lias i- 



which leads to the conclusion that zenith dry range re- 
fraction can hi* determined simply hy measuring the 
surface pressure. 

B Standard Atmosphere 

\n isothermal atmosphere is a rather p«»or assumption, 
however; a much better idea would Id* to assume a stan- 
dard atmosphere profile, For instance, a typical profile 
(see Ref. 2. pp. 82-88 1 for a station at a height h. alxive 
sea level would In a constant lapse rate atmosphere to 
II km. defined as follows: 


<)<*<(!! h) 


T(Z) - T. - yZ 


P(Z) 



y lapse rati*. K/km 


T(z) T. 



T, 21fl.fi5 K 
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and .in iMitlicnn.il atmosphere from II km to 21 km. <l< 
lined a* follows: 

(11 *.)£«£ (21 *..) 


T(») T, 



I Minj tin- above atmospheric model. the calculation nl 
dn /enith range refraction proceeds as follows; 



which is. ol course, ul< ntical to the result obtained for 
the previously considered isothennal atmosphere. 

C. Hydrostatic Equation and the Perfect Gas Law 

Imally, it (.m Im ipnte easily shown that the above 
result can easily In obt lied by merely assuming th< 
hydrostatic equation and the prrfixt gas law; 

ill' (hydrostatic equation) 

where; 

fi density 
and 

/’ i>HT (perfect gas law) 

< hie then has; 



so that 


Tlie upper limit of the second integral (21 /i„) is allowed 

to go to infinity (because of the small contribution alxivc 
21 km) so that 





S'"*©' 


10 C I 



or, once again. 

the above results lead to the unmistakable conclusion that 
dry zenith range refraction is simply a linear function of 
surface pressure: 

10 

and not of surface rcfractivitv, as was fonnerly believed 
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III. Determination of Wet Zenith 
Range Refraction 

In the previous section, It was seen that for thy zenith 
range refraction one hasi 

Mir 10 ' J*ND{z)(lz 

10 ■«'.(!) 

Correspondingly In this section, the quantity of interest 
Is wet zenith range refraction! 



where: 

PJLz) = 610 Rll(z) exp ?|~ ) 

Obviously, before one can attempt the determination of 
the Integral of wet refractivity, one must know some- 
thing about the functional dependence of temperature 
nnd relative humidity upon height above the station (s). 
Subsection lil-A below will describe typical temperature 
and relative humidity (altitude) profiles, while Subsection 
III-B will present an integration of wet refractivity based 
on several simplifying assumptions. Subsections 1II-C nnd 
1II-D present expressions for wet range refrnction ns 
derived by C, C, Chao and P. S. Callahan at the Jet 
Propulsion Laboratory in 1973. Subsection III-E presents 
more recent work (1974) by Berman on wet zenith r:.,igo 
refraction, while Subsection III-F compares the expres- 
sions for wet zenith range refraction derived in Subsec- 
tions III-B through 1I1-E. 

A. Temperature and Relative Humidity Profiles 

The functional dependence of temperature and relative 
humidity upon nllltudo obeys no simple laws that would 
lead to convenient, explicit expressions for T(z) and 
R]I(z)i about the best one could expect to accomplish 
would be to find the best approximations for T(z) and 
Rll(z) that arc not so cumbersome or restrictive that they 
effectively preclude their usage for the purposes intended 


here. A brief description of temperature and relative 
humidity profiles ensues. 

1, Temperature profiles. The region of interest In de- 
fining the functional dependence of temperature upon 
altitude Is from the station height to the tropopnusc, or 
roughly 0 to 11 km, At the tropopause, the temperature 
is approximately ' 55° C, and at this temperature wet 
refractivity is nil. It is fortunate that many of the Deep 
Space Stations (DSF ) me meteorologically similar, i,c„ 
temperate, semiarld, and not in proximity to any large 
bodies of water. For climates such ns these, one can 
Identify two basic types of temperature profiles: a 
summer-type profile (warm and usually clear), nnd a 
winter-type profile (cold nnd often cloudy), The summer- 
type profile Is the simpler of the two, nnd is shown In 
Fig. 1. 

♦ 

Generally, over a period of days, the temperature lapse 
rate and 7’, (“extrapolated surface temperature”) remain 
reasonably constant, although near die ground there is 
a diurnal swing from night inversion to day ground 
heating, anti back, The winter-typo profiles are more 
complicated in that Inversions occur up to a far greater 
altitude, and rather than one main segment, they consist 
of several segments, ns seen in Fig. 2, 

To formulate a general expression for T(z), which 
would account for the local surface effects, would be an 
almost impossible task; however, if one could extrapolate 
the temperature lapse rate down to the surface nnd define 
an "extrapolated surface temperature (T t ),“ one could 



Fig. 1. Summer typo profile 
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T~ 


Fig. 2. Wlntor typo profile 

easily construct n simple genera] function of temperature 
up to the* tropopnuso, ns follow, 


t 


Fig. 3, Rotative humidity profiles 

Since no patterns are really disceninble, one might ad- 
vnntngeously assume: 

IUI(z) ~ constant 

and since the greatest contribution to wet range refrac- 
tion occurs near the surface, where; 



T t extrapolated surface temperature 


/i //((,) ~ lill. 


T(z) Tr - yz 

.. ( r s - r, i 

v(lM i ii - k, ( 


and 


l\II, ■ surface relative humidity 


This approximation should work reasonably well for 
summer-type '.rofilos, but less well for winter-type pro- 
fllcs; however, the summer wet range refraction values 
are frequently 5 or more times the size of winter wet 
rnngo refraction values, so at least the approximation 
would be most accurate when range refraction is at its 
largest. 


2. Relative humidity profiles. After examining a num- 
ber of relative humidity profiles, one is forced to conclude 
that there appears to exist no particular functional depen- 
dence of relative humidity upon height above the station 
(such as exists with temperature)} perhaps the most that 
can ho inferred is that relative humidity seems (but only 
tenuously!) to remain more constant with altitude during 
summer-type weather than during winter-type weather — 
although this would once ngnin have a good implication 
in the far larger magnitude of summer range refraction 
values as compared to winter range refraction values. 
Some typical examples of relative humidity profiles are 
shown in Fig, 3. 


it would seem appropriate to allow; 

RU[z) - nu. 

3, Actual temperature and rclntive humidity profiles, 
Appendix A presents 10 actuni temperature and rclntive 
humidity profiles measured at Edwards Air Force Base 
during 10GS-19G9. The profiles are alternate day-night 
cases, and were chosen during the following months; 

(1) December 

(2) February 

(3) April 

(4) August 

(5) September 

to provide seasonal variation. These cases will be utilised 
to test the various hypotheses advanced in the pursuit of 
a wet range refraction model, as detailed in the sub- 
sections that follow. 
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B. Direct Integration ot Wet Zenith Refractivity 

In Subsection A, expressions for T(z) and MI{z) were 
postulated ns follows! 


where! 


T(sr)»T,-y* 

T t ~ extrapolated surface temperature 

v=3 \ir znz) 


and 


MI(z) * MI, 

Assuming wet refractivity is nil by the tropopausc, the 
integral one wishes to evaluate is: 


A/l,„ = 10* l C,C 3 


and sincei 


P»(*) 


P,„(z) = 610 R/7(a) cxp(^|^) 


then: 


A R,„ = 10-»CxC 8 


RH(z) /AT(z)~B\ . 

TO? exp \ T(*)-c )‘ L 


where 

Now since 
Rff(ss) = MI, 

then 

A R w = lO' 1 CiCjKJ7 a 


C 2 = 810 C a = 2934100 


[T(; 


1 /AT(s) - B\ , 

(»)]* exp l, T(z) - C / 


Finally, utilizing 

T(s) = 7V - 
dT -• — y dz 


ono 


hns 


A/l t o « lO-'C.Cu/l//, / ~ jexp (^~^)} [ “ “”] 
Jr, " ' r ' y 

f 10-'C,C 5 Jl//a f"l ( ... /AT-«\K„„ 

= 1 7 } 


thus one wishes to evaluate the integral 

’r, 


,„ mS f' ( AT — B\ (IT 
W(T Z ) = / cxp^-Y^-j-j^r 

Jt , 


IF the substitution 


Y 


AC — 6 
T — C 


is made, one then hns 

/AT — B\ ,, , vv 
exp = ax P + 1 ) 

= exp (A) exp (Y) 


and 


T = 


AC - B + CY 


dT = 5 y f C dY 


d? _B- AC [" Y s 1 
T- ~ Y* rfl L( AC - B + Cy ) = - 


If one defines 


Y 0 


AC-B 


then 


dT 

T 


r If dY n 
= “ B - AC L (1 + Y/Y 0 )" J 
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( 


) 

i 


so t)mt 


An asymptotic series solution is postulated ns follows; 


r f#r li\ (I T 

'V(- T 77E-)W 

J Tt 


• virii 


exp (A) l exp (V) 

n - ac / (i » r/y„y 

./ 1 ir#i 


■ilY 


For the most extreme possible range of temperatures, one 
iias for the variation in the Y-dependent terms 

216 31 6 

—22.7 SYz5 -14,5 
1.4 X 10*" , Sexp(Y)SS X 10 7 
1.15 S{l + V/y„)23 1.22 

since the variation in the denominator is almost nil com- 
pared to the variation in the exponential term, one might 
that an approximate solution could bo of the form 


/ 


exp (V) 
v7vT7 d) ■ 


exp (Y) 


(i + y/y«y 


(i -f y/y„r 


This motivates an attempt to determine a solution by first 
assuming that 


f exp (Y) 

J (i + y/v„r 


(IY 




4- Y/Y„)» 


and then attempting to solve the resultant differential 
equation in F(Y). One begins by differentiating both sides; 


ev.p ( i) 


exp (V) 


(1 Y/y„) J {1 + y/Y„) 


-f(V) 


(±\ esp(Y) m 

\yJ (l + y/Y 0 ) 3 ( } 

exp (Y) dF(Y) 

(1 4- Y/y„) a clY 


so that 


1 •- F(Y) - (JLY — 
(> \yj (i + y/y n 


dF(Y) 
) clY 


or 


mr) 

clY 




f<Y) - Y — 4 — . 
'* 4<i + Y/Y„)< 


t " 
« 


(IF(Y) vx A) ( \( \ \ 

~HT 4 ( 1 4 y/y„)/*i y >)\ x) 

i 0 

_ ( i \ A; 

4\V../{1 + Y/Y») f " 


To deduce a recursion relationship, the terms in F(Y) are 
manipulated: 


u/v\ . . V A j 


F(Y) = >" 




(1 4- Y/YnV 


a,. 4 y 


^ (i 4- y/y,)' 

i i 

; A +V ^ 

" W (1 4 Y/y„) ( *' 


and 


(y„) (i 4 y/y„) F{Y) 

I y_A_ 

\yJ (i + y/y») 4 (i + y/y 0 v 

*J±\V 4 

Vy.Jwfi + y/y,,)^ 


so that 

CC 

_v4 


V ( i \ A; 

W \Y„/ (1 4 * Y/Y 0 )W u 

« M 

V Aj n / _2_\ y Aj 

+ 4 (i 4 - y/Y c )M u’o/4 


j n 

which requires 
Ao = l 


(i + y/y„)^ 


-1 = 0 
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mill 


Aj 


V) / 

V ) A/u 

4 i(i + y/y„ ) M ‘ \v„; (i + 

> U 1 


diverges ns / -» <*>, viz: 

(/ + !)»" term (/ + 2)!(Y 0 )'(1 + Y/Y 0 )' 


Aj 


Y„/(i + r/y»v* , i 


so that 


unci 




i + 2 


W v ~ : ’ 

Y 0 (l + Y/Y„) 


(/)•" term (Y 0 )'+»(l + Y /Y 0 )' + ‘(/ + 1)1 

/ + 2 

Y u (l + Y/Y„) 


> 1 when | / + 2 1 > | Yo + Y | 


or 


Ij.i — A; 


j’ 2 + / \ 

l Y„ 


/ 


Combining A 0 = 1 unci the above recursive relationship, 
one lias 

A„ 1 



with the final result that 

pm = V <Ldi4>! 

4 (YoV(l + Y/YoV 


and 


««P(Y) f/y exp(Y) va (f + l)l 

v/v \2 /i v/v ft (V 


(1 + Y/Y 0 )* (1 + Y/Yo)« (Y 0 )J(1 + Y/Y 0 )' 


Now the series 


o (/ + 1 ) ! 


24 


/cO 


(Y.)'(l + Y/Y 0 )' 


However, a solution is still valid for this asymptotic series 
as it converges very rapidly for small /: 

/ = 0 term = 1 

i = 1,or,, ’ = WTT 
3 -2 


f c: 0 tOlTll ■■ 1 ■ — ■ - pp -4- 0,000 

I - teim f y/ y u) , ~ l 


so that 


/ 


eX P (Y) ,y 

V/V M! J 


2 ! 


31 


(1 + Y/Y.)“ 

- - - 0X Pi Y J_. I] 4 , 

' (1 + Y/Y,,)- Y 0 (l + Y/Y„) (Y„)»(l + Y/Y„)« 

+ ..] 


or 


J (L + Y/Y„)» 


exp (Y) 

(1 + Y/Y,,)* 


The final derived expression for wet zenith range refrac- 
tion then becomes 


AH,, 


^'C.C. M . 1 W(T>) 


( w'C,c,nn,\ T'l i ( at - b Yi r ,, r 

= { 7 1 1 F {“'(—)} 

( 10-*C l C,nr/.exp(A) 

1 y(B“A 


( 


«p(Y) , 

1C) f | (1 + Y/Y 0 )- 




10-' C.C.R/i, exp (A)) 

exp (Y) 

y(B - AC) J 

1 + Y/Y,,)'- 


r<r,i 
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and since 


TSk statistics of Table 1 above are ns follows: 


exp [TO] 

'exp [TO] 


3 X 10-' 


one can ignore the upper limit, so that 


A/!„i 


{ 


10-'C,C s /?/7,exp(A)l 
y(l! - AC) j 


exp [TO] 

[1 + TO/**]* 


♦■rmuhi — 5*5 cm 
ir^ny — 1,4 cm 


•rioui = 4,0 cm 


bins,i min = +4.4 cm 


Now 


binsdny “ — 0.G cm 


exp (A) exp [V(T f )] 



and 


[i + to/W - } J 


with the final result that 


. „ ao-'CiC./u/, ] r/, cyi / at ,- b \ 
A/ " { y<« - AC) / Lv ~ Tj ) J exp ( T, ~ c ) 


The ten cases listed in Appendix A were numerically 
integrated (from this point on to be defined as the 
“actual’') and compared to values computed from the 
nbove equation (now to be defined as the "Bemrinn 70”), 
The results are seen in Table 1, 


Although this model works reasonably well for the day 
cases, It is obvious that the combination of strong night- 
time temperature inversion and ground distortion of (an 
increase in) relative humidity causes the model to be 
ineffective for night enscs, Considering the composite 
<r — 4,0 cm, one would be better off using monthly aver- 
ages (and with far less trouble!) of AR U . than one would 
be using this model, according to the work of V. J. 
Ondrasik and K. L. Thulecn on usage of monthly wet 
refraction averages (see Refs, 3 and 4), As a sidelight, 
however, S. C. Wu (Ref, S) Ims been successful in apply- 
ing the model above 3,000 m, which is quite reasonable 
in view of the fact that it is the near surface distortions 
in temperature and relative humidity that make the model 
ineffective for surface measurements, 

C. The Chao Model 


, In 1973, C, C, Chao produced the following model 

Table 1. Berman 70 vs Actual (“Clino”) to predict wet zenith range refraction (see 


Ref, G): 


Case No, 

Herman 70, 
cm 

Actual, 

cm 

A, cm 

J J 1 a 

Night cases 




1 

13.3 

4.8 

4-8.5 


3 

11.9 

3.8 

4-8.1 

Chao began by assuming the following: 

5 

5,3 

3.7 

4-1.0 

~ [PYK] = — [piu]g; hydrostatic equation 

7 

21.0 

18.1 

4-2.9 

0 

Day cases 

10.8 

9.7 

+1.1 

PW = pto[fl\Y / ]T; perfect gas law 

2 

3.7 

4.0 

—0.9 

T a* T„ ~ yzj temperature lapse 

-I 

4.5 

2.9 

+ 1.0 


a 

3,5 

4.0 

-1.1 

where 

8 

17.0 

19.3 

-2.3 

pto — water vapor density 

10 

5.4 

5.7 

-0.3 

R\V = perfect gas constant for water vapor 
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The nbovo three equations led to unreasonable results, 
however, and Chao, considering the perfect gas Jaw in- 
applicable, replaced It with the ndiii 1 tic laws 

PW - W[pw]P 

where 

f) = specific heat ratio (. — '1,3 for water vapor) 

In actual fact, there is no reason to question the adequacy 
of the perfect gas law for water vapor so long as conden- 
sation does not occur, As P. S. Callahan quite correctly 
pointed out (sec Ref. 7), the real problem lies with a mir- 
npplicntion of the hydrostatic equation. For a mixture if 
gases, Dalton’s Law (see Ref, 2, p. 18) states that the 
total pressure P is equal to the sum of the partial pres- 
sures of each constituent gas (/); 

r*=£ej 

so that for the atmosphere, the hydrostatic equation Is 



with 

p* = weighted density 

but it is not necessary that for each individual constituent 
gas 

particularly in the case of water vapor, which (relatively 
speaking) has a very low saturation vapor pressure and 
is immensely affected by local surface effects, i.e., bodies 
of water, etc. Although the derivation of the model was 
flawed, the use of a constant parameter to fit actual data 
gave the model validity, and, in fact, it represented quite 
an improvement over the only other surface wet model 
then existent (Berman 70), Table 2 presents a comparison 
of the Chao model with the data for the ten test cases in 
Appendix A, 


Table 2. Chao vs Actual 


Caso No. 

Chao, cm 

Actum, 

cm 

A, cm 

Nljjht eases 

1 

2,4 

•1.8 

**2,4 

3 

2.7 

3.8 

- 1.1 

5 

2.0 

3.7 

-1.7 

7 

1-1.7 

18.1 

a.-t 

0 

3.3 

0.7 

0.4 

Day cases 

O 

1.5 

■1.0 

“3,1 

4 

2,-1 

2.0 

-0,5 

0 

3,0 

<1.0 

-1.0 

8 

21.9 

19,3 

+ 2.0 

10 

5.1 

5.7 

-0,0 


The .statistics from Table 2 arc as follows; 

"■in Bin ~ 3.5 cm 
— 1.9 cm 
mniui = 2,8 cm 
biasninht ~ “3,0 cm 
bias,i liy ~ —0.5 cm 

These can be compared to Chao’s own published figures 
(from Ref, 6): 

o’iiibIii — ‘hi cm 
(r,| ny - 2.0 cm 

The performance of the model in regard to day cases 
is creditable, although even the day residuals appear to 
have a definite negative bias; for instance, in Fig. 4, p, 41, 
of Ref. 0, 17 of 19 day case residuals are negative. For 
night cases, however, the model does poorly; even more 
troublesome than the relatively large standard deviation 
is the veiy large negative bias (—3 cm). This very large 
night case bias is likewise readily seen in the same Fig. 4, 
Ref. 6. The overall performance of the Chao model will 
bo compared to other models in Subsection III-F, 
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where 


■j 1/1 + 0.078 T, (°G)i fit parameter 


■ DAY NW, 

WITH "AVERAGE ' 1 
PROFILE 


-jj- = 8.7 + 0/13 T, (°C){ fit parameter 

U ~ height at which water vapor vanishes 
o r* 

erf (*) « exp ( - {/■) (hi 


xXX 

s. \ 


NIGHT NW, 
WITH "AVERAGE" 
PROFILE 


The basis of this model was an empirical fit to water 
vapor pressure (data from western Europe) in the form of 

FW(z) ~ I’W, exp {-as - bz?) 

in addition to the use of a constant lapse rate temperature; 

T(z) = T, ~ y z 

which, when combined, simply yields 


Fig, 4. Day and night "average" wet refractivlty profiles 


D. The Callahan Model 

In late 1973, P. S. Callahan derived a model to predict 
wet zenith range refraction, ns follows (see Ref. 7)i 




/ P\V,cxn(-<iz-hz-) , 

Art,, - lO-C.C J Uz 


The model in its full form is rather awkward; Callahan 
indicates that for the following nominal values; 

a * 0.218 knr 1 
b « 0.018 knr 1 


y — 7 K/km 
II 10 km 


T = 300 K 


“ Crf (iV^)] '' 7Wf[ CXP (“ 4b) 


The model reduces to ('‘Callahan”) 1 


. „ _ us x io-*pw, 

A iw ~ (T,/300)- 


This model, now to bo called “Callahan,” is compared to 
the Appendix A test cases in Table 3. 




'll Is to be understood that the model loses accuracy for 7'„ outside 
the range 200 K < T, < 310 K. 
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Toblo 3. Cullnhnn vs Actual 


Case No, 

Cullalmn, 

cm 

Actual, 

cm 

A, cm 

Night cases 

1 

3.8 

4.8 

-1.0 

3 

4,1 

3,8 

4 0.3 

5 

3.2 

3,7 

"0.5 

7 

15.4 

18,1 

-2.7 

0 

4.8 

0.7 

-4.0 

Day cases 

o 

2.0 

4,0 

-2.0 

•i 

3.7 

2.0 

4-0.8 

0 

5.0 

4.0 

40,4 

8 

21.2 

10.3 

4-1.0 

10 

0.5 

5.7 

+0,8 


The statistics of Table 3 are ns follows: 
o»u:,i = 2,6 cm 
ir,i„ y = 1.3 cm 
n-totm - 2.0 cm 
bias,, id, i = — 1.8 cm 
biaSdnj. = 4> 0.4 cm 

The ir,i n j value can be compared to Callahan’s own pub- 
lished day case standard deviation (from Table 2, “All," 
Ref. 7): 

(j.i’iy = 1.4 cm 

The above results are substantially better than the Chao 
model in both the standard deviation and the bias; addi- 
tional analysis of the Callahan model will be presented in 
Subsection 1II-F. 

E. Empirical Approach to the Prediction of Wet Zenith 
Range Refraction 

In the course of developing a new radio frequency 
angular tropospheric refraction model (Refs. 8, 9, and 10), 
this author found it necessary to determine the ratio of 
wet zenith range refraction to dry zenith range refraction, 
i.e., to determine some f such that 



At that time it was empirically found that there existed 
strong correlation between the ratios of wet and diy 
zenith range refraction and of wet and dry surface ref roc- 
tivily, i,e„ 



with 

K~0,32 

There was an immediate implication that, if indeed, there 
existed strong correlation, this (assumed) relationship 
might be useful in predicting wet zenith range refraction. 



in 1970, J. V. Ondrnsik considered a very similar 
approach, that is, correlating zenith wet range refraction 
with surface wet refroctivity (Ref. 3, p. 34): 

AR, 0 = K„ [NW t ] 4* Kj 

but (apparently) dropped the idea as unpromising. 

The 10 cases in Appendix A were uliliz-d to perform 
a least squares curve fit to 

Afl w vsK[|~].W. 

which resulted in the following value of K: 

K = 0.3224 
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or, lilt 1 following model (“Berman 74"); 


Tablo 4. Borman 74 vs Actual 



In Table <1, this model is compared to the Appendix A 
test cases. 


The statistics of Table <1 are as follows! 

u’nurtit ™ 2,6 cm 
ir.iny - 1.6 cm 
” 2,2 cm 

biaSntRiH ~ “1.9 cm 
bias, i nr “ + 0,6 cm 


Case No. 

Herman 7 A, cm 

Actual, cm 

A, cm 

Night cases 

t 

3.5 

4.8 

- 1.3 

3 

3.8 

3.8 

0.0 

5 

3.1 

3.7 

-0.0 

i 

15.5 

18.1 

~2.0 

0 

•1,0 

0.7 

-5.1 

Day enser 

o 

2.5 

4.0 

“2.1 

4 

3.0 

2.0 

+ 0,7 

0 

5.1 

4,0 

+ 0,5 

8 

21.0 

10.3 

+2.0 

10 

0.8 

5,7 

+ 1.1 


As can be seen from a comparison with Table 3, the 
above model and the Callahan model produce substan- 
tially similar results (about which more will be said in 
the next Subsection). 

Since all the models (with the exception of Herman 70) 
show u very strong negative bias for night cases when 
compared to day cases, one is motivated to consider the 
causes of and possible solutions to this difficulty, The 
main considerations that detract from an “average” wet 
refrnctivity profile arc fluctuations in relative humidity; 


profiles. This in fact turns out to be the case. Once again, 
a model with the following form is postulated; 


A R w = K(RH,) 



where 

[K,i day profiles 

K ~ l 

l K , , night profiles 


Rii(z)^nii, 

and near surface variations in temperature. As lias already 
been commented on, the fluctuations in relative humidity 
are difficult to categorize, while the near surface tem- 
perature undergoes a fairly regular diurnal swing, ns can 
be seen in both Figs, 1 and 2, 

If one assumes an “average” wet refrnctivity profile, 
then obviously using a “night” N\V, will give too small 
a total wet zenith range refraction; conversely, using a 
“day” N\V t will lead to too large a total wet zenith refrac- 
tion, ns is (schematically) seen in Fig. 4. 

The conclusion one expects is that the K( — 0,3224) 
determined for a mixture of day and night cases would 
be larger for night only cases and smaller for day only 


A least squares curve fit was performed on the Appendix A 
test cases. This process yielded 

K, - 0.2896 
K, ~ 0,3773 

or 

1 2519 [i—i] exp (4 y^ct); night profiles 

This model (“Berman (D/N)”) is compared to the Appen- 
dix A test cases in Table 5, 


exp gj : day profiles 
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Tablo 5. Borman (D/N) vs Actual 


Citso No. 

Berman (D/N), 
cm 

Actual, cm 

A, cm 

Night cases 

1 

*1.1 

*1.8 

-0,7 

3 

4.4 

3,8 

+0.0 

5 

3,0 

3,7 

— 0.1 

7 

18,1 

18.1 

0.0 

0 

5,4 

0,7 

-4,3 

Day cases 

0 

2.2 

4,0 

—2.4 

4 

3.2 

2,0 

+0.3 

0 

4,0 

4.0 

0,0 

8 

10,7 

19.3 

+0.4 

10 

0.1 

5,7 

+0.4 


The statistics of Table 5 arc ns follows: 
nniRiit = 2.0 cni 
o*du j. ^ l.X cm 
uiotni — 1*0 cm 
biaSniR,,, = -0.9 cm 
bias,i„j. = —0,3 cm 

As can be seen by comparing the above with Table 4, 
the separate day-night model represents n definite im- 
provement over the composite model, 

As an alternate attempt to explore methods to account 
for the systematic diurnal surface temperature variations, 
the following was tried. 

Let 

T m m ss lowest previous 24-h temperature 
T max ~ highest previous 24-h temperature 

Then, to moderate the night and day temperature profile 
distortions, define: 

T (night cases) =* 


T (day cases) eE-Sl+Isli 

A new K for the above defined model (now to be called 
the "Herman (TMOD}") wns calculated via least squares 
ns: 

K ~ 0,3281 

This model is compared to the Appendix A test cases in 
Table 0, 


Tablo G, Borman (TMOD) vb Actual 


Case No, 

Berman (TMOD), 
cm 

Actual, cm 

A, cm 

Night cases 

1 

*1.0 

4,8 

-0,2 

3 

4.8 

3.8 

+ 1.0 

5 

4.1 

3,7 

+0,4 

7 

17,5 

18.1 

-0,0 

0 

0,3 

9.7 

-3,4 

Day eases 

o 

2.0 

4.0 

-2.0 

*i 

2.9 

2.0 

0.0 

0 

4,0 

4.0 

— 0.0 

8 

19,8 

19,3 

+0.5 

10 

5.3 

5.7 

-0.4 


The satisfies of Table 6 are ns follows: 

Unlfilu = 1*0 cm 
cr,|ny — 1.2 Cm 
Wtotnl = 1*4 cm 
biasniRhi = —0.6 cm 
bias, i a y = —0,6 cm 

As can be seen by comparison with Table 5, the results 
are exceedingly similar to the separate day-night 
model, as one would expect because of the similarity in 
assumptions. 

Finally, Subsection III-F will compare the various 
models discussed in this Section. 
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F. Comparison of Wet Zenith Range Refraction Models 

Tliy results of each of the models previously presented 
nre summarized in Table 7, 


Table 7. Composite model comparison 


Parameter 

Herman 

70 

Clmo 

Callahan 

Herman 

74 

Herman 

(D/N) 

Herman 
('I'M OD) 

°nl|!lii 

6.5 

.3,5 

2.0 

2,0 

2.0 

1.0 


1,-1 

1.0 

J.3 

1.0 

1.1 

1.2 


4,0 

2,8 

2.0 

O <1 

1.0 

1,-1 

l)lns„j ( , tll 

*4.4 

—3,0 

-1.8 

1.0 

-0.0 

-0,0 


0.0 

H),5 

4 0.4 

4 0.0 

■0,3 

-0.0 


Additionally, the model residuals* for the test cases are 
seen in Figs. 5 and 0, while Fig. 7 compares the Chao, 
Callahan, and Berman 74 model for RII, = 0,3 and a 
temperature range of 270 K to 310 K, 


22 
20 

\L 
16 
14 

6 12 
% 

< 10 
8 
6 
4 
2 
0 

Fig, 5. Test case comparison of the Berman 74, 
Callahan, and Chao models 
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8 
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• R f 


ta 


■ 2 

A □ 


(....I .1 I I I . 1 . 1 I 


234567B9I0 
CASE No. 

Fig, 6, Test case comparison of tho Berman 74, 
Berman (D/N), and Berman (TMOD) models 


The models can be characterized as follows: 

(1) Berman 70: day cases work reasonably well, but 
night cases are totally useless, primarily because of 
relatively high near surface relative humidity and 
strong nighttime temperature inversions, both of 
which combine to cause the model to yield exces- 
sively high values, 

(2) Chao: day case experience is reasonably good, but 
night eases show a high standard deviation. Even 
more troublesome, the entire model has a negative 
bias, small for the day cases, but very large for the 
night cases (this is readily apparent by examining 

F ig. 7). 

(3) Callahan: day case experience is very good and 
night cases work reasonably well. Night cases still 
show a substantial negative bias, however, 

(4) Berman 74. the results of this model are very simi- 
lar to those obtained with the Callnlm.i models. 
That these two models arc very similar can he seen 
in Fig, 7; however, this similarity is not surprising 
if the models nre put in the same form: 
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f 

i 


i 
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Fig. 7. Comparison of tho Borman 74, 
Callahan, and Chao models 


Callnhan: 


und Berman 74; 




1,15 X10' 2 PW, 
(7’/300)- 


1.18 X lO'-TW, 

” (77300) 

The* difference (in cm) between (T/300)-' and 
(T/300)’ 2 only becomes significant at very high sur- 
face temperatures (T„S 305 K), 


(5) Berman (D/N): very good experience is obtained 
with both day and night cases. This model is similar 
in form to both the Callnhan and Berman 74 mod- 
els, except that it has separate coefficients for 
night versus day cases to make allowances for the 
qualitatively known effects due to diurnal tempera- 
ture fluctuations. In lief. 6, Chao presented a large 
data base of AR W versus PW, measured at local 
noon and at three locations: 

(1) Madrid, Spain 

(2) Yucca Flats, Nevada 

(3) Waggn, Australia 


By assuming fixed relative humidity lew ’ me can 
make a rough comparison of the various models 
here described, Appendix B displays this data for 
the fo flowing models ? 

(1) Chao 

(2) Berman (D/N); day coefficient 

(3) Callahan 

and for the following relative humidity levels; 

(1) 16,0% 

(2) 22,5% 

(3) 30.0% 

(4) 37.5% 

(5) 45.0% 

As an example, for the ]ll l sa 30.0% case (Fig. B-3 
of Appendix B), the models yielded the following 
statistics when passed through the data: 


Parameter 

Berman 

(D/N) 

Chao 

Callahan 

CTday. CHI 

2,03 

2.22 

2.30 

bins, cm 

+0.13 

— 0.33 

+0,99 


The above results tend to qualitatively substantiate 
the model descriptions previously presented in re- 
gard to these models. 

(6) Berman (TMOD); very similar assumption to the 
Berman (D/N) model nnd very similar results In all 
respects. 


At this time, tho best choice for an overall model, when 
considering standard deviation, bins, nnd complexity, 
would appear to be the Berman (D/N); 

1934 exp ^ * day profiles 

2519 ex P (4 ^~ c r ) ’ ni S llt: profiles 

Finally, it should be pointed out that all the models 
described here (with the exception of the Chno model) 
can be put in a very similar form, i.c., if 

AH* = PW, + ^ygyoj + py3Q0j7 
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T 

f 
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onu would Imvc for each of the following models: 

(1) Herman 70 (with T t for T,)i 
A, « 1,33 X 10* 3 

A, t= —0.3d X 10- 3 
Aj a +0.02 X 10* 3 

(2) Callahan: 

An “ 0 

A, — 0 

At «= US X 10*» 


(2) Night usago: 


The integrated dr>’ refrncttvlty is (close to) exact, while 
the wet model accuracy is considered to be: 

■ 1.5 — 2.0 cm 


(3) Herman 74: 
Ao«0 

A, a 1,18 X 10- 3 
A„ = 0 


Possible improvements to the above model (wet portion 
only) are us follows: 

(1) Use of the Callahan functional form (K[NW,]) 
with corrections for diurnal temperature effects, 


(4) Herman (D/N); 

A 0 = 0 

A, - 1.00 X 10- 3 (day); 1,38 X 1(H (night) 

A a — 0 

(5) “Chao Integrated Average Hefractlvity Profile" (see 
Refs, 0, 10, 11, and 12) a : 

A r) « 0 

A, "0 

A s - 1.08 X 10- 3 

(0) "Winn, et. nl” (see Ref. 13) 3 : 

A 0 = 0.80 X 10* 3 (footnote 3) 

Ai = 0 

A s — 0 


Since the most successful (wet) models have been 
strictly empirical, it would rather straightforwardly 
seem that a lim.'r proportionality to rcfraclivily 
would be tile uwst appropriate, i.e., consider: 


AR W *= 10* f U N\V(z)d* 


~ 10 -'NW. , 


N\V(s) 

NW, 


where 


NW(0) _ , 
NW, ' 1 


IV. Summary 

At the current timo the most accurate zenith range 
refraction model would be composed of the integrated 
dry refrnctivity (from Section I) and the Herman (D/N) 
wet zenith range refraction mod*!, as follows: 

(1) Day usage: 


*Theso two previously -....iienlloiied expressions have been added 
for the snkc of comparison. 

3 As a side light, tills coeillcleiit Is about 30% smaller than more 
frequently used duta from (Inland) (racking sites. It is conjectured 
that tills Is a result of taking surface measurements at limes when 
a near surface, high relative humidity marine layer was existent. 


NW(H) _ 

NW, 

Since one, in essence, has no other information, 
then the quantity 



fit parameter 


simply becomes tlie fit parameter, loading to the 
Callahan form. However, it is also conceivable that 
something more (empirically) representative of the 
data might bo obtained by allowing more fit pa- 
rameters in terms of tlie dominant variable, T,, in 
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which disc it might be useful to attempt a fit In the 
form: 

* n " r i' Wl { Ao + (7V300) + (r,/3W)»} 
where A 0l Ai, and A» are fit parameters. 

(2) Correction for systematic diurnal, seasonal and local 
(site) effects, 

The separate day and night coefficient model 
(Herman (D/N)) Is a step In this direction; ulti- 
mately one could visualize fitting to the following 
pnntr Hers; 

(1) Tout time of day (diurnal) 

(2) T dot ~ time of year (seasonal) 


so that a complete model could bo represented as: 


aR„ cr P\V, 


V' A|/(7'(|MT/f , W>v) 

1 * fj 

L/ o 


Even with the nbove model, this nulhor considers 
that the Information content Inherent in surface 
measurements would not allow the standard devia- 
tion of a surface meteorological measurements 
model to improve to any better than: 


ir — 1.0 — 1.5 cm f 


ns compared to the standard deviation of tho best 
current models of: 


(3) Tracking station location (I) 


tr— ' 1,5 — 2.0 cm 


Definition of Symbols 


A 

7.4475 In (10) = 17,1485 

/ 

Index variable 

A/ 

series coefficients 

K, K, 

constants 

a 

coefficient, empirical wate* vapor pressure function 

l 

length 

D 

2034,28 In (10) - 4684,1 

m 

total refrnctivity 

«i 

"Inverse scale height," km 

N. 

m 

b 

coefficient, empirical water vapor pressure function 

ND(z) 

“dry" refrnctivity 

C 

38.45 

ND, 

ND{ 0) 

c, 

0.776 


"wet” refrnctivity 

C a 

810 C 3 = 2934100 

N\V, 

N\V( 0) 

C a 

4810 

n(z) 

index of refraction 

c 

speed of light 

m 

total atmospheric pressure, N/m 3 

V 

intermediate variable 

p, 

P(0) 

f 

intermediate variable 

■P\V(s) 

water vapor pressure, N/m ! 

a 

gravitational acceleration 

P\V, 

P\V( 0) 

u 

height at which water vapor disappears 

R 

perfect gas constant 

h a 

station height nbove sea level 

RW 

perfect gas constant, water vapor 
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relative humidity (100 'r ■ 1.0) 

1 

time 

nil, 

mm 

0 

velocity 

ah 

range refraction, cm 

W 

intermediate variable 

Ml, 

"dry” range refraction 

y 

variable (AC- B)/(T « C) 

AIL 

"wet" range refraction 

r„ 

constant [AC • B)/G 

■ML 

range refraction at elevation 0 

Z 

height above station 

r 

range, km 

p 

specific heal ratio (1.3 for water vapor) 

P 

refracted range, km 

y 

lapse rate, K/km 

m 

temperature, K 

0 

elevation angle, deg 

T, 

210.05 K 

P 

density 

T, 

n 0) 

pw 

density water vapor 

r. 

extrapolated surface temperature 

p * 

"weighted'' density 

T mt 

maximum 24-hour temperature 


ft P \ T| 1 

T 

1 min 

minimum 24-hour temperature 

V 

standard deviation f 1 1/S 1 £ M I ^ J 
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Appendix A 
Test Cases 
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Fig A- 1 . Temperature and relative humidity vs altitude; Edwards AFB. 
December 9. 1968. 2 a.m. 
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Fig. A-2. Temperature and relative humidity vs altitude; Edwards AFB, 
December 9. 1968. 1 p.m. local 
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Fig. A-3. Temperature and relative humidity vs altitude; Edwards AFB, 
February 3. 1969. 2 a.m. local 
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Fig. A 4. Temperature and relative humidity vs altitude; Edwards AFB. 
February 3. 1969, 1 p m. local 
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Fig. A-7. Temperature and relative humidity vs altitude; Edwards AFB. 
August 9. 1969. 1 a.m. local 
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Fig. A 10 Temperature and relative humidity vs altitude: Edwards AFB. 
September 25. 1968. 12 a m local 
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Appendix B 

Model Comparisons to aR w vs PW s Data 
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Fig. B-2. Companion of the Berman (D/N) Callahan, and Chao models for RH 
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